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Abstract 

We study the realization in a model of graphene of the phenomenon 
whereby the tendency of gauge-field mediated interactions to break 
chiral symmetry spontaneously is greatly enhanced in an external 
magnetic field. We prove that, in the weak coupling limit, and where 
the electron-electron interaction satisfies certain mild conditions, the 
ground state of charge neutral graphene in an external magnetic field 
is a quantum Hall ferromagnet which spontaneously breaks the emer- 
gent U(4) symmetry to U(2)XU(2). We argue that, due to a resid- 
ual CP symmetry, the quantum Hall ferromagnet order parameter is 
given exactly by the leading order in perturbation theory. On the 
other hand, the chiral condensate which is the order parameter for 
chiral symmetry breaking generically obtains contributions at all or- 
ders. We compute the leading correction to the chiral condensate. We 
argue that the ensuing fermion spectrum resembles that of massive 
fermions with a vanishing U(4)-valued chemical potential. We discuss 
the realization of parity and charge conjugation symmetries and argue 
that, in the context of our model, the charge neutral quantum Hall 
state in graphene is a bulk insulator, with vanishing longitudinal con- 
ductivity due to a charge gap and Hall conductivity vanishing due to 
a residual discrete particle-hole symmetry. 
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1 Introduction and Summary 



Graphene is a 2- dimensional honeycomb array of carbon atoms whose elec- 
tronic properties at energies below a few electron volts are modeled by emer- 
gent relativistic Dirac fermions with U{4) chiral symmetry [T]-[in]- Though 
all indications to date are that clean, suspended graphene is indeed a gapless 
semi-metal, the idea that Coulomb and other interactions could potentially 
gap the electron spectrum by breaking the U{A) symmetry has inspired con- 
siderable discussion [ll]-[19] and some numerical work [20]- [29]. Moreover, 
recent experiments show plenty of evidence [30|-[l2] that some form of sym- 
metry breaking does occur in graphene at sufficiently low temperatures and 
in strong magnetic fields. Understanding the mechanism for this magnetic 
field driven symmetry breaking has become an important problem. 

It has been known for a long time that an external magnetic field enhances 
the tendency of interactions to break chiral symmetry and generate a mass 
gap in 2+1 dimensional quantum field theories with massless fermions [15] . 
This idea, which is called "magnetic catalysis", has already been applied 
to graphene and graphene-like systems [S]- [17], [11] [12]. The mechanism of 
magnetic catalysis is simple to understand [12] [38]. A magnetic field flattens 
the electron energy bands by quenching the kinetic energy to form Landau 
levels. In a state where a given Landau level must be partially filled, one 
must therefore adjust the way in which the levels are filled to minimize the 
potential energy, as all configurations have approximately the same kinetic 
energy. Then, even a very weak interaction which favors an asymmetric 
population of the single particle states in the band will break symmetry. 

An alternative, but closely related scenario for breaking the U{4) sym- 
metry of a Landau level is called quantum Hall ferromagnetism [5U]. This 
idea is inherited from multi-layer quantum Hall systems made from con- 
ventional semiconductors where it explains the appearance of ferromagnetic 
Hall states [51] [52]. The quantum Hall ferromagnet is designed to mini- 
mize the Coulomb energy of a single Landau level by making the electronic 
wave-function as anti-symmetric in the electron positions as possible, since 
this gives the electrons as much spatial separation as possible and mini- 
mizes the repulsive Coulomb potential energy. For the lowest Landau level 
of graphenci, this is accomplished by making the f/(4) indices of electrons 
occupying the Landau level as symmetric as possible which for some com- 

^We shall refer to the Landau level which apprears at the apex of the graphene Dirac 
cones as the "lowest Landau level". It is also called the "charge neutral point" Landau 
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mensurate fillings; quarter-filling, half-filling or three-quarters-filling of the 
level is a f/(4) ferromagnetic state that breaks the U{4) symmetry sponta- 
neously. If mixing of Landau levels is ignored and the dynamics restricted 
to the charge neutral point Landau level, the quantum Hall Ferromagnet is 
an exact eigenstate of the charge density operator and it can minimize the 
energy due to a density-density interaction like the Coulomb interaction. 

Quantum Hall ferromagnetism is sometimes viewed as a competitor of 
magnetic catalysis. The main distinguishing feature of the two are the order 
parameters which for quantum Hall ferromagnetism is a Hermitian matrix- 
valued U(4) electron density 

p'^'^' = (^[<n^),^^(^)]) (1) 

whereas, for magnetic catalysis, it is an electron mass operator condensat^ 

^ ^^;'ix),^;ix)]'j (2) 

In the above equations, is the second-quantized electron field operator 

(see ( !T5|) . ( !T6|) . ( IT8|) below). The repeated indices 77 in ([1]) and ([2]) are summed 
over the pseudo-spin components, 77 = 1,2 and a and a' are U{A) (spin and 
valley) indices with values a, a' = 1, ...,4. Under a f/(4) transformation im- 
plemented with 4x4 unitary matrix U, both p'^'^ and S°"°" are Hermitian 
matrices which transform by conjugation, p — )■ UpW , S — UT,W . The 
Dirac conjugate is defined as = '^'^'^'-f^ and we will use a convention (in 
(fT7|) below) for Dirac matrices where 7°^, = diag(l,— 1). The commutator 
ordering in ([T]) subtracts a constant vacuum charge density. It is impor- 
tant for giving the charge density operator the correct charge conjugation 
(particle-hole) symmetry. The constant that is subtracted is the contribu- 
tion to the charge density of the positive Carbon ions in graphene which is 
such that the particle-hole symmetric point in the electron spectrum is also 
charge neutral. This normal ordering is also helpful for defining the opera- 
tor in ([2]), as the individual terms in its computation are divergent and the 
commutator cancels some of the most severely divergent parts. 



level. 

^For a discussion of the most general form of fcrmion bilinear condensates in graphene, 
see Ref. [57] and also [SH]-[ST]. 
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However, it was observed in Ref. [12], and we will also make this point in 
the following, that in a f/ (4) symmetry breaking phase in a magnetic field, 
the order parameters ([1]) and ([2]) can both be non-zero. In this Paper, we will 
study this issue using Hamiltonian techniques and perturbation theory. We 
will concentrate on the case of charge neutral graphene in a magnetic field 
and at zero temperature. In neutral graphene, when the U{A) symmetry is 
unbroken, the four-fold degenerate Landau level which resides at the charge 
neutral Dirac point should be half-filled. It is the fact that this level is 
half-filled which is responsible for the anomalous integer Hall effect that was 
observed in early experiments [3] [9] - the charge neutral point is half way up 
a Hall conductivity step, rather than being on a Hall plateau, as it would be 
in a conventional integer quantum Hall system. 

Symmetry breaking would have this four-fold degenerate Landau level 
splitting into levels of lower degeneracy. The relevant splitting results in 
either two 2-fold degenerate levels or four nondegenerate levels, half of which 
have energy lower than the Dirac point and the other half have energy above 
the Dirac point. In the charge neutral case, the spectrum is then gapped. 
At zero temperature the levels below the Dirac point are completely filled 
and those above are completely empty. The result is a new Hall plateau at 
the charge neutral point. This plateau has been observed in clean suspended 
graphene at low temperatures and high magnetic fields [30]- [37]. As we shall 
argue, discrete symmetries imply that the bulk Hall conductivity of this 
state should vanish, i.e., the charge neutral Hall plateau is a quantum Hall 
insulator with bulk conductivities vanishing, axy = and cr^x = 0. Of course 
currents can be carried by edge states which could exist. It is also known 
that the occurrence of edge states and their detailed properties depend on 
the U{4) orientation of the condensate which we shall not address here. 

Our study will be limited to the case of weak interactions of the specific 
charge density-charge density form given in the Hamiltonian (fT5|) below, an 
important example of which is the Coulomb interaction. We will ignore 
impurities as well as explicit U{4) symmetry breaking, which would be small 
but present in graphene due to lattice effects and Zeeman and spin-orbit 
coupling. For simplicity, we will concentrate on the charge neutral Hall 
state and focus on bulk properties of graphene. We will not discuss edge 
currents which are generally important for the phenomenology of quantum 
Hall states. We shall hereafter refer to "graphene" as the U(4) symmetric 
continuum quantum field theory with Hamiltonian f|T5|) . 

The weak coupling limit has the advantage that we can definitely identify 
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the correct ground state. Our results will show that, if interactions are 
repulsive and are weak enough, in the context of our model, charge neutral 
graphene in a magnetic field indeed experiences spontaneous breaking of 
the f/(4) symmetry to U{2) x U{2). Furthermore, a particular form of the 
quantum Hall ferromagnet state proposed in Ref. [50], with perturbative 
corrections from Landau-level mixing, is the preferred ground state. The 
condensate contains a non-zero expectation value of a density and a mass 
operator, which in the leading order have identical magnitudes, but which 
have significant differences when interactions are taken into account. 

In the leading order, the condensates p'^'^ and E"^*^ have equal magni- 
tudes. This is a result of the fact that, in a magnetic field, the single-particle 
wave-function of an electron in the Landau level which resides at the Dirac 
point, i.e. the zero mode Dirac spinor, has only one or the other pseudospin 
components nonzero, not both. This is the upshot of the index theorem ap- 
plied to the continuum single particle Dirac Hamiltonian [Ijjfl This means 
that zero modes influence the matrix 



(where ri,( are pseudo-spin labels), only in one diagonal component, either 
the ri-(=l-l component or the 2-2 component, for definiteness, assume it is 
the 1-1 component. (Off-diagonal components in the indices rj, ( would vi- 
olate rotation symmetry and must therefore vanish in rotationally invariant 
states. We assume that the condensate preserves translation invariance, in 
that p°"'^ and S'^'^ are independent of position. We are also assuming that 
the U{1) subgroup of f/(4) corresponding to overall phase symmetry is not 
broken. With these assumptions, Eq. ([3]) then gives a complete description 
of the possible bilinear (dimension 2) electron condensates. If only the 1-1 
component can have an expectation value, the mass and density must be 
tuned so that p'^'^ = . We shall find that this is the case to the leading 
(zeroth) order in an expansion in the interaction strength. Higher orders in 
interactions de-tune these condensates. We shall argue that the U(4) charge 
is given exactly at the zeroth order, that is, that it is not corrected in per- 
turbation theory. On the other hand, the chiral condensate gets corrections, 

■^It also implies that a given zero mode wave-function alway has support on only one 
of the two triangular sublattices of the hexagonal graphene lattice. The zero modes in the 
two valleys live on opposite sublattices. 








(3) 
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generically at all orders in perturbation theory. Furthermore, in the rela- 
tivistic theory, these corrections can be large as, depending on the nature of 
the interaction, they can be ultraviolet divergent. 

Charge neutral graphene has important discrete symmetries, a combina- 
tion of which survives in the presence of a magnetic field and restricts the 
Hall conductivity. Given the conventions for Dirac matrices in Eq. f[T5]) be- 
low, a particle-hole transformation, C, (or charge conjugation in field theory) 
is implemented by the replacement 



C ■ 



{x) 



1 

1 



[x] 



mi* 



(4) 



This transformation maps the Dirac equation with external magnetic field B 
onto the same equation where B is replaced by —B. In addition, it transforms 
the condensates as 



(5) 



Note that both p""^ and Tj"" go to their transpose. 

Parity, P, is the reflection of one coordinate (xi, X2) — )■ (— xi, X2) and the 
Dirac spinor transforms as 



P : 



'm'l{Xi,X2) 


— ^ 


' 1" 




_^^(Xi,X2)_ 


-1 





'^?(-Xl,X2) 
^l{-Xi,X2) 



^l{-Xi,X2) 
-^^(-Xi,X2) 



(6) 



It also maps the Dirac equation with a constant magnetic field B onto the 
same equation where B is replaced by —B. As well, the condensates trans- 
form as 



(7) 



We remark that, since they interchange the upper and lower components of 
the spinor, C and P both interchange the two triangular sublattices of the 
honeycomb graphene lattice. We also remark that parity as we have defined 
it here is unorthodox in that it does not interchange the K and K' points of 
graphene, as the transformation defined in the more conventional way would 
do. Ours and the conventional one differ by a U{A) transform (see Ref. |49j 
for a detailed discussion of the conventional parity, charge conjugation and 
time reversal symmetries). 

Since the magnetic field changes sign, neither charge conjugation C in 
(jlj) nor parity P in (jH]) are good symmetries when the external magnetic 
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field is applied. However, we still expect that a combination of the two 
transformations, CP, is a good symmetry since the combination preserves 
the sign of the magnetic field. Intuitively, in a given external field, the 
cyclotron orbits of particles and holes would have opposite orientations, so 
even if particles and holes have the same energies, particle-hole symmetry is 
lost. However, we could recover the symmetry by augmenting it with a re- 
orientation of the space which reverses the directions of cyclotron orbits and 
therefore maps a hole cyclotron orbit onto a particle cyclotron orbit. This 
reorientation is parity. The combined transformation can be a symmetry of a 
system in a magnetic field. Indeed, under the transformation CP : B ^ B 
and the Dirac equation with magnetic field is mapped onto itself. However, 
the condensates that we are interested in are not automatically invariant 
under CP, instead they transform as 

CP : p""' ^ -p"'" , S""' ^ -S"'" (8) 

and having a non-zero value for either condensate would signal the potential 
spontaneous breaking of CP symmetry. We will be interested in the situation 
where these condensates are not zero, but where they also do not break the 
CP symmetry. We remind and stress to the reader that this applies to the 
charge neutral state of graphene, which we refer to as the charge neutral (or 
which is commonly referred to as the z/ = 0) Hall state only. It does not 
apply to the other Hall states since they occur at a non-zero charge density 
and non-zero charge density would break C and CP symmetry explicitly. 

To emphasize the importance of CP invariance, we note that a non- 
zero Hall conductivity would break CP and must vanish if CP is unbroken. 
One way to see this is to note that CP forbids the appearance of a Chern- 
Simons term in the effective action for the electromagnetic field which could 
be obtained by integrating out the electron degrees of freedom |53]-|55], a 
well-defined procedure once their spectrum is gapped. Absence of a Chern- 
Simons term means that the bulk Hall conductivity of the charge neutral 
Hall plateau should be zero. It is indeed observed to be zero in experiments 
[30] [M] [35] where the i/ = conductivity plateau observed in strong magnetic 
fields has axy = 0. 

If we require that CP is not spontaneously broken, it places a strong 
restriction on the form of the condensates, p^'^ and S*^"" . One possibility 
which is clear from is that they could simply be antisymmetric matrices 
so that CP : p""' -p""'" = p""' , E""' -S'^''^ = J:""' is indeed CP 
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symmetric. In fact, due to f/(4) symmetry, this is the most general constraint 
on a CP-symmetric condensate: there exists a basis where both p and S are 
Hermitian antisymmetric matrices. 

We could also use a ?7(4) transformation to choose a different basis 
where p"" is diagonal, rather than being anti- symmetric. The eigenvalues of 
an anti-symmetric Hermitian matrix are real and come in positive-negative 
pairs, so the most general form of diagonal p is p"" = diag (pi, p2, —pi, —P2)- 
Generally, E'^'^ might not commute with p'^"" and therefore, it might not be 
diagonalizable at the same time. However, in the basis where p is diago- 
nal, E'^'^ is a unitary transform of an imaginary anti-symmetric 4x4 matrix, 
which is specified by six real parameters. The upshot of this discussion is that 
the possible CP-invariant condensates are described by 8 real parameters, 
which can be taken to be the two independent eigenvalues of p and the six 
entries in the upper triangle of S*^"" in the basis where it is anti-symmetric. 

When p'^'^ has the form p'^'^ = diag(pi, p2, — pi, — P2), if pi 7^ P2, this 
condensate would by itself break the U{4) symmetry to the subgroup which 
commutes with p, the maximal torus, U{1)^. In the special case where 
Pi = P2 = P, and p°"°" = diag(p, p, — p, — p) the residual symmetry is greater, 
the symmetry breaking pattern being t/(4) — )• U{2) xU{2). This is the largest 
residual symmetry that is compatible with a non-zero condensate and CP 
symmetry. A non-zero S'^'^ could break some or all of the remaining sym- 
metry. What we will consider in the following is the scenario where this does 
not happeE0 and the system has the maximal residual symmetry U{2) x U{2) 
(we will later demonstrate that it is indeed the relevant one), i.e. p°"°" has 
the form 

p'^'^' =diag(p,p,-p,-p) (9) 
and then, to preserve the same symmetries, T,'^'^ must be of the form 

S""' = diag(S, S, -S, -S) (10) 

In summary, given the symmetry breaking pattern, U{4) — t- U{2) x U{2) 
preserving CP, the fermion bilinear condensate is completely characterized 
by the two parameters p and S in ([2]) and ffTUl) . respectively. The tuning 
that we have discussed above (following Eq. ([3])) gives one additional piece 
of information, that in the very weak coupling limit, the amplitudes of the 

^We note here that a simple generahzation of the Vafa-Witten theorem [56], which 
normally applies to vector-coupled relativistic gauge theories, to the present model of 
graphene with a Coulomb interaction implies that CP breaking would not occur. 
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two condensates are almost equal, and S ~ ±p and therefore Yf"^' ~ ip'^'^'. 
In fact, this is just what we shall find for the quantum Hall ferromagnet state 
at the lowest order of perturbation theory. For the quantum Hall ferromagnet 
state, 

1^1 / X 

S = ^ll + A,+ ... ] (12) 

We shall argue that the U(4) charge density is exact to all orders in per- 
turbation theory. The correction to the chiral condensate, Ai, is positive 
and is of linear order in the interaction potential Vix — y). The remainder, 
represented by three dots is at least of order two in the interaction potential. 

The leading terms in the above two condensates do not depend on the 
interaction strength, but only on the magnitude of the constant magnetic 
field. Furthermore, the 5-dependence must be as shown to give p and S the 
correct scaling dimensions of inverse length squared. The dependence of the 
condensates on the sign of B in ffTTj) and f|T2|) is compatible with C and P 
transformations. In the basis where they are antisymmetric, C maps p"'^ 
and S""*^ to their transpose and p'^"" obtains a minus sign. In the basis where 
p'^'^ and are diagonal, C must therefore contain a U(4) transformation 
which interchanges the positive and negative eigenvalues of both p°"'^ and 
E°"'^ . This flips the sign of diagonal p'^^ and Tj"" . This transformation, 
together with the sign flip of B under C makes the diagonal p"" and 
odd and even under C, respectively, as they should be. P maps B — )■ —B. 
This leaves p in ffTTl) unchanged and flips the sign of S in f|T2l) . as expected. 

We compute he term Ai in f|T2|) in Section 4. It is of linear order the 
interaction and it is entirely the result of Landau level mixing. Furthermore, 
it is positive and, depending on the interaction, it can contain an ultraviolet 
divergence which could make it anomalously large. This result suggests that 
the chiral condensate is much larger than the charge condensate. In fact, if 
we attempt to obtain the condensates from a phenomenological description 
of fermions where we make the approximation that they are non-interacting 
and we add terms to Hq with a U(4) chemical potential J (Px \1/'^^Pq°" 
and mass term j (Px ^'^Mq" where, 

pI"' = diag {po,po,-po,-po) , Mq"' = diag (mo, mo, -mo, -mo) 
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we would find that the best fit has /iq 
theory are then 





\B\ . , 


Po = 


— sign(mo 








—sign (mo) 



1 + E 



0. The condensates in the free 

(13) 
(14) 



which is very similar to ( ITTj) and (fT2|) . Note that, in (fT4|) . the second term in 
Eo is formally ultraviolet divergent and it is entirely due to the contribution 
of higher Landau levels. It would thus be reasonable to interpret our 
result for the condensates in the interacting theory to be entirely due to 
mass generation. The mass that is generated breaks the U(4) symmetry 
to U(2)xU(2) and it preserves C and P invariance (and the magnetic field 
further reduces C and P to CP). 

The model of graphene that we shall use has four species of relativistic 
electrons with f/(4) symmetry and an instantaneous interaction 



E 



If 

H H ^;U^)V^c(^) + 2 / ^'"^^'^ p{x)V{x-y)m^) 

^ 1 ^ A 1 ^ 



(7=1 r),C=l 



where is the non-interacting Hamiltonian given in the first term and 
Hi is the interaction Hamiltonian in the second term, '^^^[x) is the second 
quantized electron field obeying the anti-commutation relation 



vi/;;(x),vi/ft(^-')|=5-'5^^5(f_^-') 



The single-electron Hamiltonian has the form 



h 











Here, the conventions for relativistic Dirac matrices are 
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1 
-1 



7 



1 
-1 



, 7 







(16) 



(17) 



(18) 



'^It has meaningful interpretation once the divergent part is deah with by multiphcative 
renormahzation. 
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The Hamiltonian fll7p describes electron dynamics at energies near the 
Dirac points [T]. We shall use a convention where both graphene valleys 
have the same Hamiltonian. This convention assigns the upper and lower 
components of the spinors in each valley to the two graphene sublattices 
differently - the upper/lower component in one valley live on sublattice A/B 
and the upper/lower component in the other valley reside on sublattices 
B/A. In this way, both valley and both spin states have the same single- 
particle Hamiltonian and the emergent U{A) symmetry is explicit. We will 
work in units where ^ = 1 and the emergent speed of light, the graphene 
Fermi velocity, vp = ^- The covariant derivative, Di = di — iAi, couples the 
electrons to a constant external magnetic field B, and we could use either 
the symmetric gauge A{x) = ^ {—y, x) or the Landau gauge A{x) = B (0, x), 
for example. The charge density is given by the trace over U(4) indices of 
the operator in ([1]), 

p(^) = EE^ (19) 

(T=l 77=1 

The interaction potential V{x — y) is assumed to be instantaneous and trans- 
lation invariant. We shall also assume that it is a positive definite kernel, 
that is, that its Fourier transform V{k) = J (Pxe~'^^'^V{x) is non-negative, 
V{k) > 0. This means that the potential is repulsive at all distance scales. 

Here, we have used the continuum, low energy theory of graphene which 
has explicit U{4) symmetry. Of course this symmetry is not exact but is 
emergent at wavelengths much larger than the lattice scale. In the action 
which describes our model, the free field term 

So= [ dt [ d^x li'^'-f^D^ij" 



has S0(2,l) relativistic invariance as well as U(4), whereas the interaction 
termqj 

r . r ,0 ^ 1 



^Coulomb = -7— / dt I d^x ^^^r^^r'^r' 



^The nonlocal nature of the Coulomb interaction is due to the fact that the electro- 
magnetic fields extend out of the graphene plane into three space dimensions where 

Vcoulomb (X) = = ^ 6 (X) (20) 

47re|a;| Aire J\j2 

The Coulomb interaction is approximately instantaneous since the speed of light in vacuum 
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is U(4) but is not relativistically invariant. In spite of this, it is known that, 
in the absence of background magnetic field, the field theory with action 
5* = 5*0 + 5'couiomb is renormalizable in the sense that ultraviolet divergences 
encountered in diagrammatic perturbation theory can be canceled by adding 
counter-terms with the same form as those operators already appearing in 
the action. The operators ii/j'^'j^Dtip'^ and itp'^^ ■ Dip" obtain different coeffi- 
cients and the speed of light has logarithmic renormalization and becomes a 
scale-dependent parameter where it runs to larger values at low energies and 
momenta [62]- [67]. The graphene fine structure constant agraphene = 4^^^;^^^ , 
which controls the strength of the interaction, then runs to smaller values in 
the infrared, making the weakly coupled field theory stable under renormal- 
ization. The Coulomb interaction in graphene is thought to be strong, posing 
a difficulty for the validity of this perturbative approach, which would by 
valid if, for example, there were some screening mechanism which makes the 
Coulomb interaction weaker than the naive coupling Ographene = ~ fif ■ 
It is fair to say that there is presently no quantitative picture of how such 
a screening mechanism would work. In this Paper, we shall also take the 
assumption that the interaction is weak and we explore the consequences. 
Though being perturbative, our techniques are not quite as elegant as the 
Feynman diagrams which are used in Refs. [62]- [67] since we must first deal 
with the vacuum degeneracy which occurs at the leading order in perturba- 
tion theory. 

The interaction in ( JTSll was modeled on this Coulomb interaction which, 
when only components of the fermion operators which have wave- vectors very 
close to the K and K' degeneracy points are kept, have the U(4) symmetric 
form with charge density given in f|T9|) . It is possible to keep track of the 
leading explicit f/(4) symmetry breaking terms in the Hamiltonian, which 
arise from the lattice scale Coulomb interaction. They appear as higher 
dimensional local operators that are suppressed by powers of the lattice con- 
stant. This was done by Alicea and Fischer [68]. It is also possible that 
there are other short ranged interactions that are compatible with the lat- 
tice symmetry and which do not have the form of a charge density- charge 
density interaction that we are considering. An important example is the 
electron-phonon interaction which has recently been discussed in Ref. [69] . 

is about 300 times larger than the Fermi velocity, which is the emergent speed of light. The 
Coulomb potential contains no dimensional parameters and is therefore scale invariant. 
The quantum field theory defined by (fT5|) with V the Coulomb potential ([20l) is scale 
invariant at the classical level. 
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Generally, in the continuum limit, such interactions are also suppressed by 
positive powers of the lattice spacing and are therefore small. We will ig- 
nore all of these corrections while we study the problem of spontaneous ?7(4) 
symmetry breaking in the model f|T5|) . The symmetry breaking terms, as 
well as some others, for example the Zeeman interaction which, due to the 
presence of the external magnetic field, breaks the spin symmetry can be- 
come important later, in that they can determine the U{4) orientation of 
the condensate once the symmetry is broken. This can, for example, favor a 
truly ferromagnetic state with all spins aligned with the magnetic field over a 
sublattice symmetry breaking antiferromagnetic state. In spite of important 
differences in their physical properties, in the t/(4) invariant model, these 
states would be related by f/(4) rotations and would therefore appear to be 
indistinguishable. 

In this Paper, we are going to find that, if the interaction strength is 
sufficiently weak, and when the interaction is positive and couples charge 
densities as in f|T5|) . the quantum Hall ferromagnet state, defined in Eq. 
and denoted by |qhf), is the ground state of ffTB]) . |qhf) is an eigenstate of the 
single-particle Hamiltonian Hq and has a particular distribution of electrons 
in the half-filled Dirac point Landau level, and with all Landau levels below 
the Dirac point filled and all above the Dirac point empty. When interactions 
are taken into account, |qhf) is not the exact ground state. It is corrected in 
perturbation theory by the usual formula, 

V If 

|gs > = |qhf > - — — - / £x(fyV{x - y) p{x)p{y) |qhf > + ... 

Hq — Z J 

(21) 

which mixes it with higher Landau levels (the second term in fl^T]) is a su- 
perposition of single-particle states, all of which have particles and holes 
occupying higher Landau levels). 

In the following, we prove that the state fl2Tl) is the ground state of f lTSj) 
when the strength of the interaction is sufficiently weak that perturbation 
theory is accurate. This proof includes the contributions of all Landau levels. 
This is an extension of previous results where this was demonstrated in a 
lowest Landau level approximation [50] [52]. Indeed, for the scale invariant 
Coulomb interaction, all Landau levels do contribute to the interaction energy 
at the same order of magnitude and Landau level mixing could be a strong 
effect. 
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2 The Quantum Hall Ferromagnet State 



In this section, we will construct the quantum Hall Ferromagnet state using 
the second quantized formalism. First we need to study the single particle 
states in a magnetic field. In the symmetric gauge, where A{x) = {y, —x), 
the single-particle Hamiltonian ( ITTll can be written as 

h = 



V2Ba^' 




(22) 



where we have assumed that B > and we have introduced oscillator oper- 
ators 



a 



IX 



a 



^[id, + dy-^{y + ix)\ (23) 
^[id,-dy + ^{y-ix)\ (24) 



with non-vanishing commutators [a, a^] = 1 and = 1. Note that (3 and 

do not appear in the Hamiltonian. They contain the guiding center coor- 
dinates. The single-particle energies are obtained by solving the eigenvalue 
problem 

hipnmix) = Elpnmix) (25) 

The oscillator ground state is the Gaussian function 



/o(^) = foix^y) 



2tx 



/3/o = 



(26) 



We use this function to construct an infinite tower of zero modes of (125 



Um{x) 















(27) 



and positive energy (electron) and negative energy (hole) Landau level states 



±^^{n-l)m(^). 



n= 1,2,...; m = 0,l,... , E = ±V2Bn 



where 



(at)" (/3t)m 
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fo{x) 



(28) 
(29) 



The second-quantized electron field is, 

oo oo oo 

^''{X, t) = J2 + E E i^^ri^i^Xm + ^^nTJixKl) (30) 

m=0 m=0 n=l 

The non-vanishing anti-commutators of zero mode, electron and hole creation 
and annihilation operators are 

(31) 

respectively. With completeness sums for wave-functions, these guarantee 
that the second quantized electron operator fl30l) obeys the anticommutator 
(ITB]) . The non- interacting part of the Hamiltonian is diagonal in the single- 
particle states. Its second quantized form is 

oo oo 4 

^0 = E E E + &nL&n™) + ^0 (32) 

n=l m=0 tT=l 

Note that the zero mode creation and annihilation operators do not appear 
in this expression. Eq denotes the ground state energy, which is equivalent 
to the sum over zero point energies of the electron and hole states. The total 
electric charge is given by 

/4 oooo4 oo4^ ,\ 

dx^^p^^i^) = E E E - KlKm)+Yl E i^^rlK - 2 ) 
cr=l n=l m=0 a=l m=0 o-=l ^ ^ 

(33) 

As usual, electrons and holes have positive and negative charges. The zero- 
modes on the other hand carry fractional charge ±| per degree of freedom. 
This is familiar from zero modes in other contexts, such as fractionally 
charged solitons [70] [71], and was discussed in the context of graphene in 
Ref. [1], but will not play an important role here, except to illustrate that 
charge neutrality requires that half of the zero-mode states must be filled. 

In perturbation theory, which we shall employ to study the full Hamil- 
tonian Hq + Hj, we begin by studying the degenerate neutral ground states 
of Hq. a ground state of Hq has all of the non-zero energy particle and 
hole states empty, that is, a ground state is annihilated by the particle and 
hole annihilation operators and 6^^. In order to be electrically neutral, 
it should also have half of the zero modes occupied. This leads to a large 
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degeneracy of the ground state corresponding to all choices of half-filling of 
the zero modes. Our task will be to study the leading order in degenerate 
perturbation theory to determine which of the half-filled states has the lowest 
energy once interactions are taken into account. 

In the next Section, we shall prove that, amongst all possibilities, a can- 
didate for a ground state that it chooses is the "quantum Hall ferromagnet" 
(qhf) state, which is one of an infinite family of states, an example of which 
is 



Here, |0) is the state with no excited electrons or holes, aj^„ |0) = 0, |0) = 
and the zero energy Landau level empty, |0) = for m = 0, 1, 2, ... and 
n = 1, 2, .... To half-fill the zeroth Landau level, we must create two electrons 
for each label (m). The creation operators in flMjl do just that. Another way 
to specify this state, which avoids the infinite product of creation operators, 
is to have it satisfy 



<^ |qhf) =0,b-^^ |qhf) = , C='''^ |qhf) = , r='^' |qhf) = (35) 



Note that we have chosen to create states with U{4) labels 1,2. In fact, we 
could have chosen any of a family of states which are related to this one by 
U{4) transformations. The most general such state is 



where, "^^^^z \caa'\'^ = 1- Here, c^a' is independent of m. For each m, there 
is an identical superposition of the six basis vectors > (with the 

six possible pairs of non-identical indices {a, a')), which are a basis for the 
six- dimensional fundamental representation of the U(4) algebra represented 
by the Young tableau which has a single column with two boxes. 

If the |qhf) state is the ground state, it breaks the f/(4) symmetry to a 
residual U{2) x U{2), those transformations that are in U{A) and which still 
leave the state in flMj) unchanged. The zero-mode parts of the generators 
of the U{2) X U{2) residual symmetry subalgebra of U{4) which annihilate 



|qhf) = n i^'X) |0) 



(34) 



m 




(36) 
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|qhf), and therefore generate unbroken symmetries are 



Q = ^l^'l + C€ + €C + C^C, Q = -^l^'l - + + €C 

(37) 

where m = 0,1,2,... is summed in each term and we have omitted similar 
terms with particle and hole creation and annihilation operators in each 
generator. 



3 Degenerate Perturbation Theory 

We have studied some of the kinematical properties of the |qhf) state. Now, 
we must establish that, if the interactions are weak enough, |qhf) is indeed 
a ground state. We shall do this by demonstrating that it saturates a lower 
bound on the expectation value of the interaction Hamiltonian. This shows 
that it is one of a family of degenerate ground states which are related to each 
other by U (4) transformations. We note here this is already well known in an 
approximation where the mixing of Landau levels is ignored [51] . The novel 
aspect of our demonstration is that it takes the mixing of higher Landau 
levels into account. The fact that the contribution of higher Landau levels 
cancels is essentially a consequence of the special symmetries of the charge 
neutral state of graphene and will not apply to the generic quantum Hall 
ferromagnet states of bilayer semiconductors. 

The Hamiltonian f|T5|) . written in terms of creation and annihilation op- 
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erators- 



H = EoiB) + J2 {ala^m + KlKm) + \ j d^x<fyp{x)V{x - y)p(f)(40) 

nma 

We will find it convenient to split the charge density operator into three 
parts, 



a=l 

where pi{x) is quadratic in 'dm and 'd]^, 

oo 4 



(41) 



pi{x} = 

mm'=0 

P2{x) is linear in 'dm and 

oo oo 4 



o-cr' = l 



(42) 



n=l mm'=0 aa' = l 



(43) 

and psi^x) is independent of ^m and t^J^, 

(44) 



P3(-) = E E E 

nn'=l mm'=0 crcr'=l 



^The zero-point energy of the non-interacting Dirac electrons can be defined by Zeta- 
function regularization, 

Eo^-j^V2C{-^)\B\iv (38) 

Here, V is the volume of the two dimensional space. We have defined the infinite sum over 
energies of Landau levels as the analytic continuation of Riemann's zeta function, 



C(-i) = lim^^n" 

■^~^^2 n— 1 



.2078862250. 



(39) 



Note that the zeta function is finite at s = — ^. The \B\^ behavior of the energy density 
is due to the scale invariance of massless fermions. Once the negative sign of the zeta 
function is taken into account, the ground state energy p8p is positive. 
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where : ... : denote normal ordering. 

Perturbation theory is comphcated by the fact that Hq has a degener- 
ate ground state. Eq is its smallest eigenvalue and the eigenvalue equation 
Hq |c) = Eo\c) is satisfied by any state |c) which is annihilated by both parti- 
cle and hole annihilation operators, aj^^ \c) = and \c) = for all m, n, a. 
These states can contain any configuration of occupations of the zero modes. 
For the following argument, we need not restrict ourselves to charge neutral 
states only. 

As is usual in degenerate perturbation theory, the degeneracy of these 
states is resolved by finding the eigenvalues of the matrix (c| Hi \c') in the 
space of degenerate ground states of Hq. Eigenvalues of the matrix (c| Hi \c') 
are the energies of the lowest eigenstates of H = Hq + Hi to the leading 
order in degenerate perturbation theory. Of course, finding the complete 
spectrum of eigenvalues is a difficult problem which we shall not be able to 
solve here. We will, however, be able to identify the ground states. By the 
standard variational argument, the eigenvectors with the lowest eigenvalues 
can be found by searching for a state |c) for which (c| Hi \c) is a minimum. 
Such a state is guaranteed to be a ground state, or one of a set of degenerate 
ground states. 

If |c) is any of the possible ground states of Hq, it is easy to see that 
(c| Hi |c) = ^ J (fx(fyV{x-y) (c| p{x)p{y) \c) 

= ]-! dVyV{x-y) J2 (c| Pi{x)pi{y) \c) (45) 

i=i 

where Pi{x) are the terms in (HTl) . given in Eqs. (H2l) . (H3l) and (jHj). We can 
treat the terms with (c| pi{x)pi{y) \c) separately for each i. 

Let us begin with {c\ pi{x)pi{y) \c) , where pi{x) is given in Eq. fH2|) . 
First of all, it is straightforward to see that |qhf) is an exact eigenstate of 
the density pi{x) with zero eigenvalue, 

pM |qhf) = (46) 

so that I J cPx(PyV{x — y) (qhf| Pi{x)pi{y) |qhf) = is indeed the minimum 
if V{x — y) is a positive kernelf]. 

® If p\{k) = ^e~'^^'^ pi{x) is the Fourier transform of pi{x), and V{x) = 
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To see that pi(x) |qhf) = 0, we observe that, when we operate pi(x) on 
the |qhf) state defined in we obtain 

Pl(f) |qhf) = Y.i^l{x)^r^,{x) (C^C' - l^mm') n 

mm' <T=1 ^ ^ m" 

(48) 

To evaluate this expression, it is easiest to separate the summation over m 
and m' into two parts, the first consisting of the terms where m = m' and 
the second consisting of the terms where m 7^ m'. 

Consider the terms with m = m'. The state |qhf > is an eigenstate of 
each term {O^^}"^ — |) with eigenvalue | for a = 1, 2 and eigenvalue — | for 
cr = 3, 4. When summed over a, these contributions cancel, for each value of 
m. Therefore, the terms with m = m' have vanishing contribution in (I48p . 

Now, consider the terms where m ^ m! which are of the form ^^^'(9^/. 
They annihilate a zero mode in state m' and create one in a different state 
m, but without changing the a-label. However, because of Fermi statistics, 
this is only possible when the m'-state is occupied and the m-state with the 
same cx-label is unoccupied. The |qhf) state is specifically constructed so that 
all states with a given cr-label are either completely filled or completely empty, 
so this process is not allowed in the |qhf) state and the terms with m ^ m' 
must therefore also vanish in (148|) . Thus we conclude that pi{x) |qhf) = 0. 

If we truncated the interaction to the lowest (charge neutral) Landau 
level, terms containing P2{x) and p^{x) would be dropped from the interaction 
Hamiltonian and we could stop here. We would conclude that |qhf) is an 
eigenstate of the interaction Hamiltonian Hj with zero eigenvalue, saturates 
a lower bound on the energy and it is therefore a ground state. We will now 
include the mixing of the non-zero Landau levels. This mixing is due to the 
presence of P2(x) and pz{S) in the interaction Hamiltonian. We will show 
that, even when this mixing is included, |qhf) remains the lowest energy 

i j cPxd?yV{x-y) {c\h{S)pM |c) = (c| j d?kp,{-k)V{k)h{k) |c) = 

^ V{k) > (47) 
e " 

where {|6)} are a complete set of intermediate states. This inequahty is saturated by 
|c) = |qhf), since, as we shah show, pi{x) |qhf) = 0. 
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E 



d^k 



(e|pi(fc)|c) 



state. 

To this end, we have aheady shown that | J (Px(fyV{x—y) (qhf | pi{x)pi{y) |qhf) = 

and saturates the lower bound on the positive semi-definite expectation 
value I / (fx(PyV{x — y) {c\ pi{x)pi{y) |c) where |c) is any of the array of 
possible ground states of Hq, that is, any state which has no particles or 
holes in non-zero energy Landau levels and has any fractional filling of the 
charge neutral point, zero energy Landau level. We will now show that 
the other two terms in the expectation value of the interaction Hamiltonian, 

1 / (PxV{x-y) (c| P2{x)p2{y) \c) and | / (fx(PyV{x-y) (c| p?.{x)h{y) |c) are 
independent of which state |c) is used to take the matrix element, as long as 
it has no excited particles or holes. This will be a result of CP symmetry. 

Consider the interaction with density P2{x) which, when we operate on 
|c), produces states with a single particle or a single hole in a higher Landau 
level, 

P2{x) \c) = i^^^\x)i^^,{x)a±r^. \c) + ^L(^)^Lna:)C'^nL \c) 

mm'n mm'n 

(49) 

and 

j d?x£y V{x - y) (c| p2{x)p2{y) \c) = 

= I SxcPy V{x -y)Y ^i.'(^)^itH^)^it^^(l/)^™'(l/) (c| J^C'C' \c) 

nmm' cr 

+ J d\Sy V{x -y)Y i^l'{x)i^l2{xWn2\y)i^m'{y) (c| J^C'C' Ic) 

nmm' a 

(50) 

Now, we note that, because 'ipl^i{x)ip'^{x) = ■j/',^, (x)'?/'im(x), the wave-function 
contributions from intermediate particles and holes are identical in the above 
equation. This allows us to combine terms so that the zero mode creation and 
annihilation operators in fl50|) combine as an anti-commutator. The result is 

j d^xd^y V{x - y) (c| p2{x)p2{y) \c) = 

= 4 J d'xd'y Vix -y)Y ^^Ux)^ljitlix)i^LtlHy)i^m'iy) (51) 

nmm' 
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The right-hand-side of this equation is independent of the state |c) which is 
used on the left-hand-side. This means that this contribution to the energy 
is the same for any |c). 

Now, consider the third contribution from P3{x) which, when operating 
on |c), creates a particle-hole pair in higher Landau levels, 

P3(x)|c)= Yl ^1^2Hx)^IjUx)<IkL'\c) (52) 

mm'nn' 

The contribution to the energy is proportional to 

£x£y V{x - y) (c| p3{x)p3{y) |c) = 
= Ajd'xd'yVix-y) ^iw(^)^itH^)^it^Hz/)^i7i.'(y) (53) 

nmn'm' 

This result is also independent of the population of zero modes in |c). 
Finally, we have demonstrated that 

(ci^Hc>4/<.w.>.,.-„x:(cip,,.)p.(.)ic>. 

= -J cPx(fy V{x — y) (c| pi{x)pi{y) \c) + { independent of |c)} (54) 

Further, we have shown that the first term on the right-hand-side is minimal 
when |c) = |qhf). This establishes that |qhf) is indeed the state with lowest 
energy in degenerate perturbation theory. 

Finally, we remark that |qhf) is an eigenstate of the Hamiltonian only at 
the zeroth order of perturbation theory. Once it has been identified as the 
state with the smallest expectation value of the energy, we use it to construct 
the ground state that is corrected to first order in perturbation theory, 

|gs > = |qhf > - — — - / Sx(fyV{x - y) p{x)p{y) |qhf > + . . . 

-no — -C/O ^ J 

which we quoted in ( l2Ti) . where V is the projection onto states orthogonal 
to |qhf > and the dots indicate contributions of quadratic and greater order 
in the interaction strength. The interaction Hamiltonian creates states with 
electrons and holes in higher Landau levels. 
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4 Condensate 



The expectation value of the U(4) charge density in the perturbatively cor- 
rected quantum Hall ferromagnet state (12T!) is 



p""' = < qhf|^ |qhf > 

- [ (fy(fzV{y- z) Re < qhf — ^— p(y)p(£)|qhf > + 

J Hq — Eq 

(55) 

and the chiral condensate is 
S'^-' = < qhf|^'"'(f)^"(f)|qhf > 

- [ (fy(fzV{y- z) Re < qhf — ^— p(jr)p(^)|qhf > + . 
J Ho — Eq 

(56) 

The leading terms on the right-hand-sides of (l55l) and (!56|) are governed 
strictly by the zero modes and can easily be computed, 



^ oo 

p'^'^' = -diag(l, 1, -1, -1) J2 = diag(l, 1, -1,-1) 



m=0 



\B\ 



1 °° B 

S--' = -diag(l, 1, -1, -1) J2 + - • • = diag(l, 1, -1, -1)4^ + - • ■ 

m=0 

which are the results quoted in ( ITTl) and ( fT2i) at the leading order. To find 
this leading term, we have used the completeness sum quoted in Eq. ( 158!) 
below. 

The second terms on the right-hand-side of (l55l) and (!56|) arise from in- 
teractions. The perturbative corrections to the condensates in fl55p and (1561) 
are quite different. We can argue that the linear order correction to the U(4) 
charge density, the second term on the left-hand-side of Eq. f l55|) . vanishes. 
This is due to the fact that, because of the projection operator, the only 
intermediate states that contribute are those with a particle-hole pair where 
both are in non-zero Landau levels, or a particle or hole in a non-zero Lan- 
dau level and a displaced zero mode. In both cases the wave-functions which 
contribute are orthogonal - they appear as ipkm\x)iljlt\x) or iljiml\x)ilJs{x), 



23 



and their conjugates, where mn and rs are the Landau level quantum num- 
bers of the hole and the particle, respectively. Because of orthogonality, the 
integral of this density over space must vanish. That is, the shift of the total 
U(4) charge, as opposed to the charge density, is zero. However, because of 
translation invariance, the matrix element must be an x-independent con- 
stant. It is thus a constant whose integral over space must vanish, it must 
therefore be zero. Thus it applies to all terms at this first order perturbative 
computation of the U(4) charge density - the entire correction to linear order 
in the interaction must vanish. 

This argument can be extended to higher orders. If the U(4) charge 
density is uniform, it can be deduced from knowledge of the total U(4) charge 
which can be gotten from taking expectation values of the four operators 

Q'' — ('^nm'^nm ^ ^nm^nm) + ( ''^rn^n ^ 2 ) ^^^^^^ 

n=l m=0 m=0 ^ 

for each value of cr. These operators commute with the free and interaction 
Hamiltonians (they are generators of the U(4) symmetry) and they therefore 
commute with all of the operators which occur in perturbative correction of 
the |qhf) state. This suggests that all perturbative corrections vanish and 
the U(4) charge density which is computed in the leading order is exact to 
all orders in perturbation theory. 

The chiral condensate fl56p . on the other hand, obtains contributions from 
the linear order. This contribution is from terms with particle-hole intermedi- 
ate states where the analogous wave-functions appear as ipLn\x)'-f^tlji.P [x) = 
ipmn\x)iljls \x). Integration of this quantity over x gives Smr^ns, rather then 
zero. This allows a contribution to the chiral condensate. It should, however 
by x-independent and we can identify parts which would be identical after 
we perform an integral. In fact, an argument similar to the one that we used 
for the U(4) charge operator suggests that we can replace the expectation 
value of the mass operator by the expectation value of 



mn ^mn 



where v is the volume of space. 

To proceed, we need the sum rule (which we derive in the Appendix 

1^1. ^2 B 



m=0 



nnm{x)ul^{y) = Me--(-^)^- MipL (58) 
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where, the summation is over guiding center labels and Ln{x) are Laguerre 
polynomials, normalized so that the first few are 



Lo = 1 



L2 = l-2x + 



Li = 1 — X 



L3 = 1 — 3x + -X 



' hx' 



and 



1 d"^ 



n\ dx"- 



With this sum rule, we find 



47r 



1 poo / n _ _ 

1 + t/ drJ—ViV2^\)e-^Y. 
Jo V ' ' p=l 



(59) 



where, as indicated, corrections are at least of order two in the interaction 
strength V. 

For the Coulomb interaction 



B 



167re ^ 22p(p!)2^ 



(60) 



For large p, the summand ~ 1/ (^A"P) and the sum is logarithmically diver- 
gent. 



5 Conclusions 

In this Paper, we have shown that the weak coupling ground state fl?Il) of 
the continuum model of graphene (fT5|) in a magnetic field is the quantum 
Hall ferromagnet corrected by higher order terms. One of our results is the 
systematic inclusion of Landau level mixing. We have shown that, because of 
the residual particle-hole symmetry of the charge neutral state. Landau level 
mixing does not shift the relative energies of the possible fractional fillings of 
the charge zero point Landau level. The quantum Hall ferromagnet therefore 
remains the lowest energy state even after Landau level mixing is taken onto 
account. 

In order to understand the electron spectrum after interactions are turned 
on, we can fit the expectation values of the bilinear fermion operators that 



25 



we have found to the same expectation values in a theory with free fermions. 
There, one would expect that, as well as the terms already present in Hq, 
the Hamiltonian which describes the fermions would acquire a term with a 
U(2)xU(2) invariant chemical potential J (i^x\l/'^^(x)/ig°" (x — (y) where 

^r'=/iodiag(l, 1,-1,-1) (61) 

For reasons which we shall explain below, this chemical potential is not 
needed and it can be set to zero. 

In addition to a chemical potential, there could be a (non-local) U(2) xU(2) 
invariant mass operator J c?^x^'^(a:)M^°' [x — y)'^'^ (y) where 

71^0^-^'= mo diag(l, 1,-1,-1) (62) 

and 

The differential operators a and are defined in Eq. (123!) . Note that they 
contain only covariant derivatives, so this mass operator is formally gauge 
covariant. The spectrum of aa^ is the positive integers, n = 1,2, ... and 

327re22-(r[x+l])2v/i ^ 

Computing the charge and chiral condensates in the charge neutral state 
with the free fermion Hamiltonian with fiQ = and mo given by f l63|) will 
result in the condensates p'^^ and T,'^'^ that we have found in the theory with 
interactions turned on. Note that the induced mass is of linear and higher 

2 

orders in the Coulomb coupling making it small in the weak coupling 
limit where the analysis is valid. 

It is interesting that, what beg U(4) ferromagnetic state with a 

U(4) ferromagnetic charge density is nominally a state with dynamical mass 
generation only. The U(4) ferromagnetic charge density in this state with 
dynamical mass generation simply results from the asymmetry of the single 
fermion spectrum with a mass and in a magnetic field - when a mass is 
present, the states which were zero modes of the massless Dirac Hamiltonian 
and which occurred at the apex of the Dirac cones are shifted and now occur 
at one of the mass thresholds, either E = +mo or E = —mo, depending 
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on the sign of B. There is the further fact that has been known for a long 
time [53] [7T] , that the spectra are otherwise symmetric - all other states come 
in pairs of electron and hole states. Thus we see that the entire "spectral 
asymmetry" of the Dirac Hamiltonian comes from precisely the same states 
which were zero modes. Furthermore, the charge of the single particle state 
is given by the spectral asymmetry [7T]. Therefore, the entire U(4) charge 
is due to these threshold modes. Further, because the interactions do not 
upset the symmetry of the remaining modes, this charge seems to survive 
perturbative corrections. 

In our detailed computations, we have used a shortcut where we deduce 
the charge density and the expectation value of the chiral condensate by first 
considering the expectation values of the relevant densities integrated over 
the volume of space. The integral results in a much simpler operator and 
the expectation values are given by simple formulae. We then deduce the 
densities by dividing by the volume factor. This is justified by the fact that 
the states that we consider are invariant under translations. The expectation 
value of charge densities should therefore be position independent constants. 
The total charge would be the integrals of these constants, which must then 
be equal to a volume factor times the constant. Particularly for the charge, 
where the integrals of corrections must vanish, we deduce that the constant 
density should vanish. One must caution the reader that there are cases 
where this proceedure does not work. For example, the summation over 
Landau level integers might not commute with doing a volume integral. The 
non-normal ordered chiral charge is one example. The expectation value of 
the operator ^^[x)'^^{x) in the Hall ferromagnet state would be 

oo oo oo 

(qhf|fri(a;)v^i(x)|qhf) = + J2J2'f'nJ i^h'^Llix) 

m=0 n=l m=0 

oo oo oo 

m=0 71=1 m=0 

If we integrate over x, due to the orthogonality of positive and negative 
energy states, the integral of the second summation on the right-hand-side 
must vanish and we would conclude that the total chiral charge density is 
given by the first term. On the other hand, of we consider the density more 
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explicitly, 

(qhf I ^\x)^\x) |qhf) = J2 i^l{x)^m{x) 

m=0 

oo oo 

+ 2 5^ 5Z [^lmix)Unm{x) - mJ„_i)„ (a;)M(„_i)„ (x) 
71=1 m=0 

and then we cancel terms which have the same lower indices in the second 
summation, we see that a term with zero mode indices does not cancel and 
we obtain 

oo 

(qhf I ^\x)^\x) |qhf) = - 5^ ^l{x)^m{^) 

m=0 

which differs from the result above by a factor of one half. Further, it agrees 
with the result that one would obtain using the commutator normal ordered 
operator that we defined in ([2]) using either route, 

1 - 1 °° B 

(qhf I - [^\x), ^\x)] |qhf) = -J2 = - 

m=Q 

This turns out to have the correct symmetry properties when combined with 
the other components computed using the same technique. 

There is another argument for the non-renormalization of the charge den- 
sity, related to the fact that an induced Chern-Simons term is also not renor- 
malized beyond one loop. The latter statement applies only when there is 
a charge gap in the spectrum (as exphcit computations show[?]). To be 
concrete, consider the Lagrangian of Coulomb interacting graphene, 

1 

L = tjji'y'^duil' V'V — ^"^'j/' 

Sire ^_V2 

which, by introducing a field Aq can be written as the equivalent theory 

If we introduce gauge fields in the central torus of U(4), A"'^' = diag{Aj^, Aj^, 
and a C and P invariant mass term m'^""' = diag(m, m, —m. — m). 
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where we have written the last term so that the Coulomb interaction is 
mediated by the U(l) component of the gauge field. 

At one-loop order, integrating the fermions out of the theory induces a 
Chern-Simons term 

Scs = l^sign(m) J d^x {A^dA^ + A^dA^ - A^dA^ - A^dA^) 

in the effective action. By a straightforward generalization of the arguments 
in Refs. [72] [73], it can be shown that this Cherm-Simons term does not 
obtain corrections beyond one loop order. Furthermore, this one-loop order 
result, when used to derive the U(4) charge density, gives precisely the result 
that we have found. 

Depending on the U(4) orientation of the magnetization, the quantum 
Hall ferromagnet can have several manifestations with significantly different 
physical properties. First of all, it can be a genuine ferromagnet with the 
spins of states in the charge neutral point Landau level aligned with the same 
orientation. Recalling that the zero modes of one of the valleys have support 
only on one of the sublattices, and the other valley on the other sublattice, the 
ferromagnet has electrons occupying both sublattices with equal amplitudes. 
This is akin to a Wigner lattice. In the same manifold of ground states is the 
anti-ferromagnet where spin up lives on one sublattice and spin down lives 
on the other sublattice. Then, there are other states where the electrons 
are paired to spin singlets, but only inhabit one valley and therefore one 
sublattice. We emphasize that there is no distinction between these states 
at the level of the U(4) invariant theory, they are all in the manifold of 
degenerate ground states that we call the "quantum Hall ferromagnet state" . 
Interactions which have been ignored in the U(4) invariant continuum theory, 
but which would certainly be present in real graphene, would be important 
for distinguishing which these states actually forms the ground state. The 
energy scale of these symmetry breaking interactions is typically of the order 
of a few percent of the Coulomb energy. 



6 Appendix: Sum rule 

The wavefunction is given by 



1 («t)n(^t)- 
Unm[x) = fo{x) (65) 

v2 ynl ym\ 
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where 



'2B 



2dr — 



/2B 



2dx — —X 
X 2 



) /o 




(66) 



and we use the complex plane coordinates x = xi + ix2, y = iji + iy2 and we 
shall assume that B is positive. Recall that and ^'^ commute with each 
other. We write the kernel as 



oo „ oo ^ 

E«-(^)<-(^) =^ E r]-T[(«')"(/5')'"/o(^)][(«')'^(/3')"^/o(z/)]* (67) 

m=0 

=2^E;^ ( 2 J [(«^)"-"^/o(-)][(«^)Wo(y)]* 

m=0 



rn=0 



(68) 



r)00-|/n\m 

^^E;^ 2 (4)"(«l)"*-"^re--(«) (69) 

m=0 ■ ■ \ / 



S 1 

'27rri! 
S 1 / 1 



1 n 



)'*(Q/t)'^*e~^^^^'''^^-*'''"?^'^ 
B 



(70) 



29^. x 



2dx — —X 
2 



2^, - 



e 4 



- (a:-2/)(a;-2/)+ f {yx-xy) 



(71) 



[2; - y]" e-f(s-y)(^-y)+f(.y^-^y) 

(72) 



S 1 



27r('n)! V2 



(0 e-^(^-^")(^-^)+^(^^-^^) [25a. - B{x - y)T [x - y]' 



B 
An 



(73) 



(n) 



27r 



(74) 
(75) 



where L„ are Laguerre polynomials. 
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